On geometry-dependent vortex stability and topological 
spin excitations on curved surfaces with cylindrical 

symmetry 

V. L. Carvalho-Santos'^, F. A. Apolonio^, N. M. Oliveira-Neto'^ 

""Instituto Federal de Educagdo, Ciencia e Tecnologia Baiano - Senhor do Bonfim 
48970-000 Senhor do Bonfim, Bahia, Brazil 
^ Departamento de Fisica, Universidade Federal de Vigosa 
36570-000 Vigosa, Minas Gerais, Brazil 
'^Departamento de Qmmica e Exatas, Universidade Estadual do Sudoeste da Bahia 
4-5206-190, Jequie, Bahia, Brazil. 



Abstract 

We study the Heisenberg Model on cylindrically symmetric curved surfaces. 
Two kinds of excitations are considered. The first is given by the isotropic 
regime, yielding the sine-Gordon equation and vr-solitons are predicted. The 
second one is given by the XY model, leading to a vortex turning around 
the surface. Helical states are also considered, however, topological argu- 
ments can not be used to ensure its stability. The energy and the anisotropy 
parameter which stabilizes the vortex state are explicitly calculated for two 
surfaces: catenoid and hyperboloid. The results show that the anisotropy 
and the vortex energy depends on the underlying geometry. 
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1. Introduction 

Geometrical and topological concepts and tools are important in many 
branches of natural sciences, particularly, in Physics. For instance, the idea 
of symmetry, which is intimately associated with geometry, is a keystone for 
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studying a number of fundamental properties of several physical systems, e.g., 
the Noether theorem asserts that there is a conserved quantity to each con- 
tinuous symmetry of the associated action. Topology, in turn, is crucial for 
classifying and for giving stability to certain excitations, such as solitons, ex- 
tending objects having finite energy, and vortices, presenting a non-vanishing 
vorticity around a given singular point or a topological obstruction. In ad- 
dition, the observed vortex-pair dissociation is the mechanism behind the 
topological phase transition Vortices and solitons have been observed 
in a number of systems, such as superconductors, superfiuids, and magnetic 
materials j2l, Isj . 

Curvature effects play an important role in the characteristics of these 
topological structures. For instance, Vitelli et al have shown that in-plane 
vortices interact not only with each other, but also with the curvature of the 
substrate Curvature is also an important factor in the magnetic systems 
behaviour, in which the interaction of the out-of-plane component of mag- 
netic vortices with curved defects must cause a chiral symmetry breaking 
in its gyrotropic motion due the thin- film roughness jsi 0]. Furthermore, 
the easy-surface Heisenberg model in magnetic spherical shells predicts a 
coupling between the localized out-of-surface component of the vortex with 
its non-localized in-surface structure, associated with the curvature of the 
underlying geometry [3] and still, the smooth and variable curvature of fer- 
romagnetic nanotorus ensures the stability of the vortex for smaller radius 
than their nanoring counterparts jsj. 

In the case of two-dimensional systems, vortices and solitons can appear 
like solutions of the continuous Heisenberg Model, which has been used to 
analyse the dynamic and static properties of vortices, showing that the en- 
ergy of these excitations is closely linked to the geometrical properties of the 
surface jof — 14|. It has also been shown that, for simply-connected surfaces, 
the vortex energy presents a divergence, which can be controlled by the in- 
sertion of a cutoff in the region where the continuous limit of Heisenberg 
Hamiltonian has not validity. In the case of magnetic systems, this diver- 
gence must be controlled by the development of an out-of-plane component 
in the vortex core region. Soliton-like solution has also been considered in 
the above cited works and it has been shown that its characteristic length 
depends on the length scale of the surface. For finite surfaces, fractional/half- 



soliton solutions have been found 11, 15 . Furthermore, the interaction of an 



external magnetic field with Heisenberg spins on a cylindrical surface yields 
a 27r soliton-like solution, inducing a deformation at the sector where the 
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spins are pointing in the opposite direction to the magnetic field 16|. A 27r 



sohton has also been predicted to appear on curved surfaces with cylindrical 
symmetry, provided the magnetic field is coupled with the curvature of the 



substrate 17, 18 



In this paper, we study the anisotropic Heisenberg Model on curved sur- 
faces with cylindrical symmetry. We are interested in studying a class of 
topological spin textures on these manifolds in such way that both, soliton 
and vortex-like solutions are considered. Solitons are predicted to appear on 
an infinite cylindrically symmetric surface, if the isotropic case is taken into 
account. In our assumptions, the soliton characteristic length is rescheduled 
to one and does not depend on the characteristic length scale of the surface. 
For finite surfaces, fractional solitons, which have not topological stability, 
are found. 

Our analysis includes the study of the XY model and vortex-like solu- 
tions are considered. It is shown that, for non-simp ly connected manifolds, 
the obstruction of the surface ensures the vortex topological stability due the 
removal of spurious divergences appearing in the core region. The XY model 
can also yield a helical-like state, however topological arguments must not be 
used to ensure the stability of this spin configuration. The energies of these 
spin textures are calculated and we we get that it depends on the surface cur- 
vature. Furthermore, we calculate the critical anisotropy parameter for which 
the vortex appears as the ground state and we show that it is also associated 
to the geometrical properties of the surface. We calculate explicitly the vor- 
tex energy and the critical anisotropy parameter for two different geometries: 
the catenoid and hyperboloid, which are negatively curved and non-simply 
connected manifolds. The choice for studying these surfaces is associated 
with the fact that both can be realized in fluid interfaces provided with an 
orientational ordered phase as a conseguence of the interplay between surface 
tension and orientational elasticity |19i]. Furthermore, catenoid is the shape 
minimizing the curvature elastic energy, appearing in phospholipid vesicles 
of high topology [20| and the symmetry of these geometries allows us to com- 
pare our results to these found for the cylindrical case, largely explored in 
nanomagnetism researches. 

To proceed with our analysis, this work is organized as follows: in Section 
|2] we present the continuous anisotropic Heisenberg Model on rotationally 
symmetric surfaces. The results and discussions for the isotropic Heisenberg 
Hamiltonian and the XY model are also considered in this section. Section 
in] brings the discussions about the model on the catenoid and hyperboloid 
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surfaces and compare our results with that obtained for the surface of a 
cyhnder. Finally, in the Section HJ we present our conclusions and prospects 
for future works. 



2. Continuum Heisenberg Model on Curved Surfaces 

The anisotropic exchange Heisenberg model, for nearest neighbour inter- 
acting spins on a two-dimensional lattice, is given by the Hamiltonian below: 

^latt = - J' J]KmJ + m\m] + (1 + A)m^m^^], (1) 

where J' denotes the coupling between neighbouring spins, and according 
to J' < or J' > 0, the Hamiltonian describes a ferro or antiferromagnetic 
system, respectively, rrii = {mf,mf,ml) is the spin operator at site i and 
the parameter A accounts for the anisotropy interaction amongst spins: for 
A > 0, spins tend to align along the internal Z axis (easy-axis regime); for 
A = 0, one gets the isotropic case; for — 1 < A < 0, we have the easy-plane 
regime, while the A = — 1 case yields to the so-called XY model, which has 



been considered on curved surfaces 2l|. If we focus on a two-component 
spin, imposing niz = 0, so that rhpRM = {fnx,my), we get the planar rotator 
model (PRM). 

In the continuum approach of spatial and spin variables, valid at suffi- 
ciently large wavelength and low temperature, the model given by Eq. ([1]) 
may be written as follows (J = J' /2): 

AJ f f 1 
H = -— / / ^=(1 + \ml)dr]idri2 

^'Jjt ilfha^il + SasX) VWldmdm, (2) 

where a is the network spacing, the surface has curvilinear coordinates r/i 
and r)2, \/M — I' 5'*"' ^^"^ are the surface and spin space 

metrics, respectively (as usual, g^^gjk = SI). Now, rh = {mx,my,mz) = 
(sin B cos sin 6 sin cos 0) is the classical spin vector field valued on an 
unity sphere (internal space), so that 6 = Q{rii,ri2) and $ = $(?7i,?72)- With 
this, the Cartesian parametrization for m yields to hat = Sab- Note that 
the first term in the first integral in the Hamiltonian ([2]) is the ground state 



4 



energy and we will renormalize it to be zero. It can also be noted that, if 
A decreases from to —1, the term depending on increases the energy if 
rriz 7^ 0, thus, the smallest energy associated to the anisotropic Heisenberg 
Hamiltonian will occur for rriz = 0. The Hamiltonian ([2]) may be also viewed 
as the anisotropic non-linear a model (NLaM), which lies on an arbitrary 
two-dimensional geometry. Thus, besides ordinary spins, the above model 
can be used to describe another condensed matter systems, e.g., a superfluid 
helium film, thin superconducting films 0, [iij, a nematic liquid crystal con- 
fined on curved surfaces [2^ or a spin lader, which consists in two or more 



coupled spin chains [24 



Our interest is to study the above model on curved surfaces with cylindri- 
cal symmetry, which, in cylindrical coordinate system, can be parametrized 
by r = {p{z), (f), z), where p{z) = p is the radius of the surface at height 
z, and accounts for the azimuthal angle. In this case, we have that the 
covariant metric elements are given by: 

1 1 

where p' = dp/dz. In this case, the Hamiltonian ([2]) can be rewritten as: 




H = J II \ [fmd.Qf + sin^ Qm?] 



+ ./^[/(e)(a,e)2 + sin2e(a,$)VW— (4) 

where /(B) = l + Asin^B. 

The Euler-Lagrange equations derived from the Hamiltonian (jlj) are eval- 
uated to give: 



2f{Q){diQ + dp) = sin 26 (9^$)' + (9^$) 



-A[(ace)V(a,ef] + -^wf-} (5) 

and 

sin^ e ((9j$ + dl^) + sin 26 {dc^Qdi^^ + 9^6(9^$) = 0, (6) 
where dC = \fgz 
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At this point, one can note that the above equations resemble, in form. 



their counterparts for the planar, cylindrical [13|, spherical [12|, pseudospher 



ical [llj and toroidal [Q"] surfaces. Indeed, whenever ( is identified with 

the Eqs. ([5]) and ([6]) recover their cylindrical, toroidal, spherical or pseu- 
dospherical analogs, respectively. Here, 6 is the polar angle in spherical 
coordinates system. In the toroidal surface, it plays a similar role, however 
6 accounts for the polar angle between the spin vector and a line parallel to 
the z-axis, crossing the center of the torus arm. TZ and s are the rotating 
and axial radius of the torus, S is the sphere radius and sinh r accounts for 
the distance measured along pseudospherical geodesic, say, a hyperbole. 

As expected, the anisotropic Heisenberg model is described by nonlin- 
ear differential equations and suitable nontrivial solutions can be obtained 
provided some conditions are imposed. Thus, special solutions for the most 
general Eq. (|5]) will be explicitly obtained by solving it for particular values 
of A. Initially, we will consider the isotropic case, which is given by A = 0. 
After, we will take A = — 1 to study the XY model. 

2.1. Isotropic regime and soliton-like solutions 

The simplest way to seek for possible soliton-like solutions associated with 
the present model on a generic surface is by considering the isotropic regime, 
A = 0, and writing down the Hamiltonian ^ and its associated equations (jS]) 
and ([H]) in a more suitable form. The assumption of cylindrical symmetry, 
9(C,0) = 0(C) and <l>(C,0) = $(0), will allow us to get the sine-Gordon 
equation in a simpler way. In this case, Eq. (jl]) is written as: 

i^isot = f d<P f ' [{d^Qf + sin2 e [d^^f] dC (8) 

J-TV JCl 

and the Eq. (E]) is simplified to: 

sin^ Qdl^ = ^ $ = Q0 + 00, (9) 

where Q G Z and 0o is a constant of integration. Other possible solutions for 
the Eq. dH]) are = and O = vr, which are the two vacua of the isotropic 
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Heisenberg model. By taking Q = 1 in the solution given by the Eq. we 
can rewrite the Eq. (jS]) as: 



2(9^6 = sin 26, (10) 



which consists in the sine-Gordon equation. By imposing the boundary con- 
ditions G(— oo) = and 0(oo) = vr, the solution for this equation is a 
topological soliton excitation, which appears as a continuous transition con- 



necting the two neighbouring minima [25 



e(C) = 2arctan(e^) . (11) 

It is important to note that, usually, the soliton presents a characteristic 
length scale (CLS) that depends on the characteristic length of the underlying 
manifold, which, in general, appears in front of the sin(2B) term in Eq. 
( !T0|) . However this characteristic length is given by g^p^/gzz, which, with our 
assumptions, is embedded in the ( parameter. Thus, the soliton would have 
its CLS rescheduled to one, e.g., if we take the cylinder case, the CLS of the 
soliton is equal to the radius po of the cylinder, however, when working on an 
infinite cylinder, the change of variable z ^ z/pq eliminates any dependence 
on Po 
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From the calculation of the energy associated to the Eq. f fTTj) . we obtain: 



E, = 2tiJ r [{d^Qf + sin^ 9] dC, 
hi 



,2C 



12) 



Ci 



which is in agreement with the saturated Bogomolnyi inequality [261 [Eg > 
SttJ), once, if we evaluate the solitonic charge, 

Q = — [ smOded^, (13) 
4vr J 

in the Eq. f fTT]) . we exactly obtain Q = — l/(e^^ -|- 1). At the interval 
( = (—00,00), the soliton agrees with its counterpart, which lies in an in- 
finite cylinder and represents a complete mapping from the spin sphere to 
the target manifold, a vr soliton, so corresponding to the first homotopy class 
of the second homotopy group of the mapping of the spin sphere. However, 
when we consider finite surfaces, such a mapping is incomplete and a frac- 
tional soliton solution appears, such that no homotopy arguments can be 
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WW 

Figure 1: [Color online] From left to right, the three first figures show the vortex pattern 
with winding number Q = 1 on the surfaces of the cylinder, hyperboloid and catenoid. The 
last figure represents a helical-like state on the catenoid. In the three first cases, the spins 
turn around the surface in a closed way. In spite of the cylinder and catenoid have different 
shape, the characteristic length associated to these are equal, then they present the same 
value for the vortex energy. On the other hand, the catenoid and hyperboloid surfaces have 
similar shapes, but they have different characteristic length and still, the mean curvature 
of the catenoid is null, while the mean curvature of the hyperboloid is nonzero. This 
leads to different values for the vortex energy on these surfaces. The helical-like state can 
present its energy equal to the vortex state, provided we have = <E>(— C). 




used for classifying solution ffTTl) as a topological excitation. Indeed, in this 
case, we must take into account the topology of the geometrical support. In 
other words, for non-simply connected surfaces, the soliton acquires a finite 
characteristic length, which prevents its collapse and, consequently, its size 
does not vanish. Similar scenarios are provided by the annulus, the truncated 



cone 
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the punctured pseudosphere [Jjj and torus [91] . 
Periodic soliton solutions can also be obtained for these rotationally sym- 
metric surfaces 13|, |27t| and the differences among the particle- like excitations 
in these and other geometries are associated to their characteristic length 



scales 27 



Now, we may wonder whether another solitonic solution, with &{(,</)) = 
0(0) and $(C,0) = '^'(C)) should not also appear in this framework. In 
fact, once ( is not, in general, a periodic parameter, sin$ would be peri- 
odic only for ( = 27rm (with m > 1 integer) and consequently, the spin 
sphere mapping should be possible only for this particular case. Besides, the 
geometry-imposed boundary conditions must lead to a periodic solution for 
6, that is, 0(0) = 0(27r), and solitonic excitations must not occur this way. 
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2.2. XY model 

Geometrically, a vortex with winding number Q may be viewed as a set 
of spins rotating in a closed circuit around a core, whose center is a singular 
point, or a topological obstruction, which makes it impossible to change the 
configuration to a perfectly aligned state without tampering with spins at 
an arbitrary distance from the core. Thus, once a vortex-like configuration 
cannot be continuously deformed to the ground state, it acquires the status 
of a topologically stable excitation. 

When we model a vortex as a continuum of spins, we intend to describe 
only its outer region; once inside the core, the analytical treatment is ex- 
pected to give only an estimate of its energy, shedding no light about its 
real structure and spins arrangement, which require numeric or simulation 
techniques. However, as we shall see in the two particular examples treated 
here, non-simply connected topologies lead to a natural cutoff for the vortex, 
which is provided by the surface hole and the solution presents no core. On 
the other hand, in the case of a connected surface, a singular core takes place, 
appearing at the surface's self-intercepting points. 

It is easy to note that a vortex solution could be obtained from the 
isotropic Heisenberg Hamiltonian, however, we have chosen to work with 
the XY model, by adopting A = — 1 in the adopted model. This choice is 
justified from the fact that it inserts constraints in the solutions for B, which 
stays confined to the x|/-plane. In this case, assuming cylindrical symmetry, 
the Eqs. ^ and ^ are simplified to: 



2cos2e9je = sin 29 




(14) 



and 



sin^ Qdl^ = 0, (15) 
and the Hamiltonian (jl]) can be rewritten as: 

HxY = J jj [{I- sin^ Q){d^Qf + sin^ Q{d^^f] dCdcj) 

+ ^ jj g"" cos^ edCdcp. (16) 

As previously noted, the simplest solution for the equations ( fT4l) and ( !T5|l 
must be given by = or B = vr. However, the last term in Hamiltonian 
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(fTB]) gives a positive contribution for the energy, and it is minimized when 
O = 7r/2, which is also a possible solution for (IT^ . Thus, in order to diminish 
the energy, the function must be confined to the xy-plane and, in this case, 
the Eq. ([15]) yields: 

$ = Q0 + 0o, Qez, (17) 

which represents a vortex with winding number Q. The winding number 
(vorticity) is formally defined, in the continuum limit, as: 

Q = ^(fm)-di, (18) 

where the integration is evaluated along a closed path C around the surface. 
The vortex pattern with Q = 1 for the surfaces of catenoid, hyperboloid, and 
cylinder can be viewed in the Fig. [1] From Eq. f|T6|) . one can see that the 
vortex energy increases with Q^, thus, vortices with Q > 1 are physically 
unstable, despite they have topological stability. Thus, we will analyse only 
the case Q = 1. 

The vortex energy, for an cylindrically symmetric surface, can be deter- 
mined from Eq. (fT6|) . and it is evaluated to give: 



(19) 



where one can note that the vortex energy depends on the ( parameter, 
which is closely related to the curvature of the surface and to the radius of 
the surface at height z. For example, the torus, sphere and pseudosphere 
present this parameter given by Ctor, Csph and Cpsph, given in the set of Eqs. 

(CD- 

Now, we can look for the possibility of the existence of helical-like states 
on the considered surfaces, which can be obtained by taking $(C, 0) = $(C)- 
In this case, the Hamiltonian yields: 

HxY2 = jj [(^ - ^i^^ Q){d^Qf + sin^ Q{d^^f] dCdcf) 

+ ^ f[ 9^' cos' QdCdcl), (20) 



and the derived Euler-Lagrange equations are: 

(21) 



2cos2e9je = sin 26 
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Using the same arguments of the previous analysis, it is easy to note that the 
simplest solution minimizing the energy is given by = 7r/2. In this way, 
the Eq. is simplified to: 

9^2$ = 0^$(C) = gC + Co, (22) 

where one can note that this solution looks like that given in Eq. (fT7|) . 
However, unlike the vortex case, q must not be an integer number, once 
( is not, in general, periodic. Thus, we have not a closed path along the 
surface to define the topological charge. In order to obtain a helical-like 
state rotating by 27i along the surface, we can assume the boundary condition 
= C)- Thus, q would be an integer and the spin vector field form a 
closed structure around the xy-plane of the internal space when the surface 
is mapped from —z to z (See Fig. [T]). The energy of this spin configuration, 
for g = 1, is given by: 

^Hei = 27rJC, (23) 

which is the same solution found for the vortex state. It is important to 
note that the helical-like state has not topological stability, once q must 
not be considered as a topological charge. The coincidence in the energy 
values of the vortex and helical-like states appears due the imposed boundary 
conditions for $(C)- Another important fact is that, if we were considering 
the magnetostatic energy term in our model, this configuration would have 
magnetic energy greater than that of the vortex state due the appearing of 
surface magnetic charges, a = m ■ n, and it would not be stable. 

Finally, once the Eq. (fT5|) gives us two possible solutions to $(0), we 
can look for the critical value for the anisotropy parameter, Ac, such that for 
A < Ac one has solution (IT7|) and for A > Ac the solution 9 = appears. 
Here, Ac will be determined by calculating the energy for these different spin 
configurations and comparing them. The stable state will be that one which 
minimizes the energy. With this assumptions, it is easy to note that 

K = -C/e, (24) 

where ^ = / ^/g^^J'g^dz and A'c = 4Ac/a^. Thus, it can be noted that, 
as well as the vortex energy, Ac also depends on the surface's geometrical 
properties and, consequently, different surfaces must present different values 
for Ac- In the case of the cylinder it is easily shown that Ac = —1. 

In order to present some examples showing the relation among curvature, 
vortex energy and critical anisotropy parameter, we will study two particular 
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surfaces: the catenoid and the hyperboloid, which are non-simply connected 
and negatively curved surfaces, having cyhndrical symmetry. 



3. Two particular cases 
3.1. Catenoid 

The first surface to be considered is the catenoid, which is a nonplanar 
minimal surface. The catenoid has mean curvature everywhere zero and has 
the fascinating property that it can be deformed into a helicoid in such a 
way that every surface along the way is a minimal surface, which is locally 
isometric to the helicoid. This surface appears in many physical systems, e.g., 
a soap film formed between two coaxial rings takes on this shape [28[ and in 
membrane fissions, where the characteristic of the neck shape is determined 



by the relationship between the neck radius and the monolayer thickness |29 

In order to describe the adopted model in the geometry of the catenoid, 
we have parametrized this surface as follows: 

X = Pq cosh I — J cos 0, y = Po '^osh | — J sin 0, (25) 

where G [0,27r], z E (—00,00) and is the radius of the circle in the 
z = plane. This parametrization yields the negative Gaussian curvature 
Ccat = ~l/Po^ech^(z/po) and it allows us to get the characteristic length 
scales of the catenoid, which are given by: 

I 2h h 
Ccat = / —dz = — , and ^cat = 2coth— , (26) 

J-h Po Po Po 

where 2h is the height of the catenoid. From the Eq. ( l26il . one can see that 
the critical value for the anisotropy parameter, in the catenoid case, depends 
on h, and it is given by: 

A'f = --cothfAy (27) 

Po \PoJ 

It can be noted that, unlike the cylinder case, the critical anisotropy 
parameter that stabilizes the vortex configuration depends on the relation 
between height and central radius (CR) of the catenoid, represented by Pg. 
For /i — )■ 0, we get A'^^* — — 1 and, for large values of h/pg, it can be noted 
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that A'|;^* decreases linearly with the height. This behaviour is expected, 
once the vortex energy, on this surface, is a linear function of h (See Eq. 
fl33|) ). The Fig. [2] show the behaviour A'^, with when we varies the height 
and CR of the catenoid surface. One can note that A'J;^* — )■ —1 for large 
values, as expected, once, when ^ h, the catenoid surface looks like a 
cylinder. 

3.2. Hyperboloid 

The hyperboloid is a quadratic surface that may be one- or two-sheeted. 
The one-sheeted hyperboloid is a surface of revolution obtained by rotating 
a hyperbola about the perpendicular bisector to the line between the foci, 
that is, about the z axis, while the two-sheeted hyperboloid is a surface of 
revolution obtained by rotating a hyperbola about the line joining the foci. 

When oriented along the z-axis, the one-sheeted circular hyperboloid with 
skirt radius p has its parametrization given by: 

X = p\/l + cos(j), y = p\/l + 1]"^ sincf) and z = hrj. (28) 

This parametrization yields a shape similar to the catenoid, however the 
last one has mean curvature null everywhere, while that one of the hyper- 



boloid is given by [30 



and the Gaussian curvature is: 

^i^yp = - 52 + (^2 + ^2)^2 • (30) 

Thus, unlike the catenoid, the hyperboloid is not a minimal surface. Fur- 
thermore, there are substantial differences in their [x, y) coordinates away 
from the plane z = 0, as well as in their geometrical properties. 

The general parametrization given in Eq. (l28|) leads to a hard integral 
to be calculated, so that, when necessary, Chyp and ^hyp will be computed 
numerically. However, a particular and interesting kind of hyperboloid, which 
leads to solvable integrals to the characteristic length parameters, is described 
by the polar hyperbolic coordinate system (biharmonic coordinates), which 
can be obtained for taking h = p m the Eq. (!28|) . This particular coordinate 
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Figure 2: [Color online] The top left figure show the critical value for the anisotropy pa- 
rameter, in function of height and CR for the catenoid, polar hyperboloid and hyperboloid 
with 6=1. Note that the catenoid and polar hyperboloid present strongly depending 
on The top right figure show for the hyperboloid with b — 0.5, 6 = 2 and 6 = 3. In 
the bottom figure, we have fixed the height z — 4 and CR varies from to 10. One can see 
that critical anisotropy parameter is highly depending on the geometrical properties of the 
surface. Furthermore, it increases with the r|a|lius, in the case of the catenoid, it decreases 
for the hyperboloidal surfaces. If we consider different heights, there was no qualitative 
changes in the A^ behaviour. 



system was recently used to develop the pseudospherical functions on an 



one-sheeted polar hyperboloid [31[. These functions may be important if 
we consider problems where the solution of the Laplace equation in systems 
with hyperbolic symmetry is demanded. For instance, they could be used to 
calculate the magnetostatic energy in magnets with hyperbolic shape or still, 
to calculate the magnetic field inside solenoids with this geometry, when it is 
traversed by an electric current and/or still, to determinate the electric field 
generated by a charged polar hyperbolic shell. 

By using the parametric equations (!28|) . the vortex energy to the polar 
hyperboloid is given by the Eq. (fT9|) with ( given by: 



Xiv) = \ , • (32) 



Cph.p = / xrnri = ^/2arcsinh(^/2,) + \ In ( | | g^^, ^ ^[I^^^^ J + 

(31) 

where 

'I + 2r/2 
l + r]^ 

and /€2 is a constant of integration. 

Finally, we must calculate the critical anisotropy parameter of this sur- 
face. In order to do this, we have gotten ^hyp numerically. The results of the 
critical anisotropy parameter for the polar hyperboloid, and the hyperboloid 
with h = 0.5, b = 1, b = 2 and 6 = 3 are summarized in the Fig. |2l As 
well as the catenoid. A' ^'^'^ depends on the height and radius of the surface. 
However, the anisotropy parameter must decreases in order to give stability 
to the vortex state. Note that for small CR, lA'J?^^! has its smallest value 
for b = 0.5, when compared to the cases where 6=1, 6 = 2 or 6 = 3. How- 
ever, when we increases the value, this situation is reversed. Finally, the 
analysis of the Fig. |2] show that A'JJ^p — t- — 1 only for — )■ and — )■ 0. 

3.3. Comparison between the vortex energy on a catenoid and hyperboloid 

From the Eq. f|T9|) and Ccat, we get the energy of a vortex on a catenoid 
surface, which can be explicitly written as: 

AnJh 

^vortex — ~ ) (33) 
Po 

which is the same energy obtained for a vortex on the surface of a cylinder 
with height 2h and radius p^. This is an interesting result, since the area of 
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the catenoid is greater than that of the cyhnder. This fact can be explained 
because these surfaces present the same characteristic length l/pg , which can 
be noted by the development of the Hamiltonian for the cylinder surface. 
Furthermore, from the analysis of the Fig. [H it can be noted that, in the z = 
plane, the neighbour spins on the catenoid must have the angle equal to that 
of the spins on the cylinder, however, when z 0, the neighbour spins that 
turning around the catenoid surface have a lower deviation one to another, 
when compared with that of the cylinder, diminishing the exchange energy 
(given by the Heisenberg model) in this plane, which compensates the largest 
area. To say this on another way, if the circle situated in the plane z = h, for 
each surface, is divided in points with a distance Ax one to another, we have 
that the circle on the cylinder is divided in n' = 27rp(,/Ax points, while the 
circle on the catenoid is divided in tt, = 2npg cosh.{h / p^^) / Ax . In this way, we 
have that n = n' cosh.{h / p^) . On the other hand, the obtained circumferences 
can be divided in arcs of angles 9 = 2ti /n and 6' = 27i/n' for the catenoid 
and cylinder, respectively. It is immediate to note that 6 = 6'/ cosh.{h/ p^). 
A closed spin texture on these circles can be given by associating one spin for 
each point, in such way that the angle between two neighbour spins for the 
catenoid (cylinder) is 6 {6'). The energy for these vortex-like structure in the 
circle on the catenoid is Ecat = = ''^cos6', where the subscripts i 

and j indicate neighbour spins. The energy for the circle on the cylinder is 
-E'cyi = n' cos 9'. Finally, Ecat and E^yi are related by: 



n cos — 



-E^cat = ir-^cyi- (34) 

n' cos ^ 



n' 



For n 1, we have that Ecat = Ecy\ cosh{h / p^) . Now, the cylinder and 
catenoid heights can be divided in m' and m circles, respectively. The total 
energy of the vortex on the surfaces are given by: 

-Etcyi = m'Ecyi and Etcat = mEcat = mEcyi cosh — . (35) 



Po 



Since m = m'^ sinh — , we obtain: 

h Po ' 



Etcat = (^^Sinhy^ Etcyl. (36) 

If — ^ 1 — 7- sinh — ^ — , and finally -Etcat = -E'tcvh as we wanted to show. 

Po Po Po 

This discrete analysis is not valid for large 2h/pg values, once in this limit, the 
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approximation of the surfaces by points does not represent the continuum 
approach to the Heisenberg Hamiltonian and the gotten results must not 
agree with that found in the Eq. fj33|l . 

In the case of the hyperboloid, analytical calculations will be done only 
for the biharmonic coordinates. In this case, in order to obtain and analyse 
the vortex energy on a hyperboloid with height 2/i, we will take the limits 
(i = —h/p^^ and C2 = h/p^ in the Hamiltonian f l2U]) . Thus, the vortex energy 
is given by: 



^vortex 



\/2arcsinh ( \f2— | H — In 



/l + 3fAV_2Aji + 2fAV\ 



^0/ V V^o 



(37) 

Despite of the polar hyperbolic coordinates give us an analytical solution 
to the vortex energy, it also will be interesting to calculate it numerically for 
the most general hyperboloid (67^ Pg) and analyse its behaviour in function 
of the radius p^. In this case, Chyp = / 'x'iji)^^^ where 



Po(i + r) 

The numerical integration was done for using the Simpsons rule with a For- 
tran code, and the results can be viewed in Fig. |3l 

To continue our analysis, it will be useful to define the upper radius (UR), 
which is the value of the radius of the surfaces in the plane z = ±/i. The 
cylinder, the catenoid and the hyperboloid have their CR and UR related by: 

=p, =pcosh(-) and i?^yP =pVlTh^, 

Upper ' Upper ' V ^ 7 Upper ' 

(39) 

where p represents the CR of the cylinder, catenoid and hyperboloid (6 = 1), 
respectively. 

From the Fig. [3], one can note that, for small values of p, the vortex on a 
hyperboloid has lower energy than that on a catenoid. However, when p in- 
creases, this behaviour changes, and the energy of a vortex on the hyperboloid 
becomes greater than that on the catenoid. The value of p for the transition 
point in which the vortex energy on the hyperboloid and on the catenoid in- 
tersect themselves depends on h, and varies from p ^ 2 for = 1.5 to p ~ 4 
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Figure 3: [Color online] Energy of the vortex on the catenoid (black line) and the hyper- 
boloid (red dashed line) surfaces. Here, we have done J = 1, 6 = 1 and = p is evaluated 
in the interval from 0.01 to 20. Four values of h are considered: h = 1.5 (top left), h — 2.5 
(top right), h — 5 (bottom left) and h = 10 (bottom right). One can see that, for small 
values of CR (central radius), the energy of a vortex on the hyperboloid is lower than 
that of the vortex on a catenoid. However, for large values of CR, the vortex energy on a 
catenoid vanishes, while that one of the hyperboloid tends to 27rln(l + h"^). The energy of 
a vortex on the polar hyperboloid (magenta dots) is always lower than that on a catenoid 
and, unlike the general hyperboloid {b ^ p), it vanishes when p ^ oo. 
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for h = 10. This fact can be explained from the differences of the properties 
of the geometries of the catenoid and hyperboloid when p — )■ oo and the 
height h is maintained constant. While the catenoid has its UR related with 
the height by cosh(/i/p), the UR of the hyperboloid increases with y/TTh^, 
which implies that, for p h, we have -Rupper ~^ -^upper -^Upper ~^ 
Indeed, for p ^ h, the catenoid surface looks like a cylinder, while the hyper- 
boloid has the topology and geometry of a plane with a hole. The analysis 
of the Fig. [3] shows also that p — > oo when -Eyortcx ~^ 27rln(l + h"^). This is 
the second time that an hyperbolic surface presents an asymptotic finite and 
nonzero energy when the radius of the surface tends to infinity [llj. This 
finite value to the vortex energy for p — t- oo can be explained because, unlike 
the catenoid case, the characteristic length of the hyperboloid does not tend 
to zero in this limit. In fact: 

x'{v),^oo = ^ =^ C£°° = HI + v'). (40) 

The case where 6 7^ 1 does not lead to qualitative changes on the vortex 
energy behaviour on a hyperboloid. In this case, the curve that characterizes 
the vortex energy is moved upwardly, however, as well as the 6=1 case, it 
tends to 27rln(l + h'^), which is associated to the fact that the b value does 
not affect the asymptotic behaviour of Chyp when p —> 00 (see Eq. (1521) ). 

Finally, we will analyse the vortex energy on the polar hyperboloid sur- 
face, which is given by the Eq. (157|1 . From the Fig. 121 one can note that 
this surface presents a vortex whose energy is always lower than that on the 
catenoid. For small CR, the vortex energy calculated by using biharmonic 
coordinates is approximately that one obtained for general hyperboloid, once 
p Ri h = 1. However, as well as the result found for the catenoid, we have 
that -Ephyp — ^ when p — )■ 00. This fact can be explained because in this 
limit, the polar hyperboloid looks like a catenoid and its characteristic length 
is given by: 

x'(^)p->oo = l^Cpi;J|r = ^/P = Ccat- (41) 

It is important to note that a vortex, understood as a topological defect, 
is a singularity in the order parameter. Then, at first sight, it could mean 
that the obtained solutions, which are smooth spin textures, are not vortices. 
However, despite we do not have analysed this case explicitly, when CR — )■ 0, 
for both surfaces, one get that the vortex energy diverges due a singularity 
in the order parameter at the center of the surfaces, which is the necessary 
condition to call these excitations of topological vortices. A way to control 
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this divergence would be to insert a cutoff of length (.^ at the center of the 
surfaces, as well as it has been done for other geometries where this singu- 
larity appeared [ol, [HI, 12|. However, unlike the spherical, pseudopherical or 



toroidal cases, if we introduce this cutoff in the point where the surfaces of 
the catenoid and hyperboloid self intersecting, these will be divided in two 
half-surfaces and this case is outside the scope of this work. 

If we consider the case where the UR of the three surfaces are equal, it 
is immediate to note that the catenoid has smaller energy than that of the 
cylinder, since the vortex energy for both surfaces is proportional to 1/p. 
However, there are not qualitative changes when we study the hyperboloid 
case. Here, as well as the previous analysis, the vortex energy on the three 
surfaces diminishes with the increasing of UR, tending to zero when r oo 
(cylinder, catenoid and polar hyperboloid). 

In conclusion, the lowest value found for the vortex energy, for any CR, 
occurs for the polar hyperboloid, which can indicate, at first sight, that 
among the geometries considered here, a "hyperboloidal nanoshell" could 
support a vortex-like magnetization with more stability than their cylindrical 
counterparts. However, to ensure this statement, one must calculate the 
magnetostatic energy for other possible magnetization states and consider 
the volume of the magnet. 

Once miniaturization is an important issue in nanotechnology, the ob- 
tained results may guide future researches on the stability of vortices in circu- 
lar nanomagnets. About this theme, it has been shown that, due the out-of- 
plane component, the vortex energy in circular nanodots is greater than that 
obtained for nanorings with same dimensions, because they do not present 



the out-of-plane component in the vortex core |32|. In this way, nanorings 
can support vortex as magnetization ground state for smaller radius than 
that of nanodots. Furthermore, it has been shown that the smooth curva- 
ture of a toroidal device must improve the vortex stability, when compared 
with cylindrical nanorings jsj. Thus, it may be possible that, by introduc- 
ing curvature in the nanorings fabrication in order to obtain "hyperboloidal 
nanorings", the vortex gets more stability, and one could fabricate smaller 
nanomagnets with a vortex as the ground state. 

4. Conclusions and Prospects 

We have studied the anisotropic Heisenberg model on curved surfaces 
with cylindrical symmetry in order to obtain a class of topological spin ex- 
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citations. By taking A = 0, soliton-like solution has been gotten and it have 
been shown that they have energy obeying the Bogomolnyi's inequahty. Frac- 
tional solitons, which has not topological stability, have also been predicted 
to appear if we consider finite surfaces. By taking A = — 1, one obtain the XY 
model, whose simplest solution is a vortex with winding number Q — 1. The 
energy and stability associated to vortices is closely linked to the geometric 
properties of the surface. Furthermore, we have gotten the anisotropy pa- 
rameter, Ac, for which the vortex appears as the groundstate, in comparison 
with the state in which the spins point in the 2;-axis direction. It has been 
shown that Ac is also related to the geometry of the substrate and it must be 
variable for surfaces with non-constant Gaussian curvature. In addition, we 
have also studied the possibility of the appearing of helical-like states on the 
considered surfaces. The energy of these states have been calculated and, for 
periodic boundary conditions, the same energy predict for the vortices has 
been obtained for this spin configuration, however, helical states does not 
present a topological charge and homotopy arguments can not be used to 
ensure its stability. 

In order to exemplify our results, two surfaces were explicitly analysed: 
the catenoid and the hyperboloid, which are non-simply connected mani- 
folds having negative and variable Gaussian curvature. For each surface, the 
vortex energy has been compared with that presented for a vortex on the 
cylinder, which is a geometry that have received large attention in nano- 
magnetism researches. It has been noted that, for small CR, the general 
hyperboloid (6 7^ p) presents a vortex energy lower than that presented by 
the vortex on a cylinder and on a catenoid, however, for p — > 00, unlike the 
other studied geometries, the vortex energy on a hyperboloid does not vanish. 
Furthermore, the smallest value for the vortex energy is given by the polar 
hyperboloid surface. Thus, it may be possible that, by introducing curvature 
in the nanorings fabrication, in order to obtain "hyperboloidal nanorings", 
vortices are more stable, and one could fabricate smaller nanomagnets with a 
vortex as the magnetization groundstate. However, when we think in nano- 
magnetism applications, this fact does not ensure that nanomagnets with 
this geometry could support a vortex magnetization configuration with more 
stability than cylindrical nanorings. To state that, one must consider the 
volume of a hyperbolic nanoring (nanomagnet limited by two polar hyper- 
boloids with internal and external radii given by pi and p2, respectively) and 
take into account the magnetostatic energy to calculate the magnetic energy 
associated to other magnetization configurations. In this way, if miniaturiza- 
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tion of magnetic elements is demanded, the understanding of the curvature 
influence on magnetic properties of nanostructures is very important, once it 
may give a way to diminish the size of nanomagnets that present a vortex as 
the magnetization ground state. In this paper, we have explicitly considered 
only the vortex state on non-simply connected and negatively curved surfaces 
by assuming that they satisfy the conditions in their dimensions ensuring this 
configuration. 

Unlike the cylinder surface, the anisotropy parameter presented by the 
catenoid and hyperboloid, which ensure the vortex as the lower energy con- 
figuration, varies with the height and radius of the surface. The results 
show that the increasing in the catenoid radius produces a decreasing in the 
modulus of the anisotropy parameter, which goes to 1 at the limit p — )■ oo, 
reflecting the fact that, in this limit, the geometry of the catenoid is similar 
to the cylinder one. The modulus of the anisotropy parameter on the hy- 
perboloid presents a different behaviour, increasing with the surface's radius 
and height. 

Despite we have given attention to magnetic systems, the adopted model 
may be relevant to study nonmagnetic systems like a superfluid helium film, 
where $ would be the phase of the superfluid order parameter 0], a super- 
conducting filrn, where the field is identified with the phase of the collective 
wave function m or nematic liquid crystals, where $ would be the nematic 



The results obtained here also may be relevant in the magnetoelastic 
membrane manipulation subject, once vortex configuration could be used 
as a way for deforming the geometry of magnetoelastic surfaces to obtain a 
particular shape. In this context, it would be an interesting task to study 
the vortex energy behaviour on surfaces with positive Gaussian curvature and 
include the surface tension term in the analysis of this problem. In addition, 
it may be important to study the influence of an external magnetic field 
on the energy calculations of magnetoelastic materials presenting a vortex 
configuration. 
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